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Abstract 



It is well known that, in the plane, the boundary of any quadrature 
domain (in the classical sense) coincides with the zero set of a poly- 
nomial. We show, by explicitly constructing some four-dimensional 
examples, that this is not always the case. This confirms, in dimen- 
sion 4, a conjecture of the second author. Our method is based on 
the Schwarz potential and involves elliptic integrals of the third kind. 
MSC 2010: 31B05, 30E20. 
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1 Introduction 

A bounded domain Q C C is a quadrature domain (in the classical sense) if it 
admits a formula expressing the area integral of any function / analytic and 
integrable in Q as a finite sum of weighted point evaluations of the function 
and its derivatives, i.e. 



where z m are distinct points in Q and a m k are constants (possibly complex) 
independent of /. 

Note: This can be generalized in various directions. One may use a 
different class of test functions, integrate with respect to a weighted density 

* Supported by NSF grant DMS-1067886. 




v 




(or over the boundary), or allow for more general distributions in the right 
hand side. "Quadrature domain in the classical sense" is used to specify the 
restricted case we consider throughout this paper. 

Suppose that Q is a bounded simply-connected domain in the plane with 
non-singular analytic boundary. Then the following are equivalent. More- 
over, there are simple formulas relating the details of each. 



(i) Q is a quadrature domain. 

(ii) The Schwarz function of dfl is meromorphic in Q. 

(iii) The conformal map from the disk to Q is rational. 



For their equivalence see [21 Ch. 14]. 

In higher dimensions, quadrature domains are defined by replacing ana- 
lytic functions with harmonic functions. We write the quadrature formula 
using multi-index notation with a = (ati, a 2 , .., at n ) consisting of nonnegative 

Q\a\ u 

integers, let I = ai + a 2 + .. + a n , and d a u = — „ nr> — - — . Then, Q is a quadra- 

<).,■;<).,■'..():;;; 

ture domain if it admits a quadrature formula for integration of harmonic 
functions u, 

p N n k 

/ udV = ^2 a m,ad a u(y: m ), (1.2) 

^ n m=l |a|=0 

where xi,X2, ..,xjy are points in Q and a m ^ a are now real constants. 

Remark: In IR 2 , any quadrature domain for harmonic functions is a 
quadrature domain for analytic functions, but not vice-versa. See [91 Ex- 
ample 1] for an example of a quadrature domain for analytic functions that 
is not a quadrature domain for harmonic functions. 

For the case of n > 3 dimensions, condition (ii) has a counterpart for- 
mulated in terms of the Schwarz potential, introduced by D. Khavinson and 
H. S. Shapiro (see Section [2]). As a consequence of Liouville's Theorem on 
the rigidity of conformal maps [3], condition (iii) does not extend to higher 
dimensions. Throughout this paper, we will make use of the reformulation of 
quadrature domains in terms of the Schwarz potential discussed in Section 

EJ 

Regarding the existence, in the case when there are no derivatives appear- 
ing in the quadrature formula (1.2), the free-boundary problem of obtaining 
a quadrature domain satisfying the prescribed quadrature formula can be 
reformulated as a so-called obstacle problem, and the existence of a solution 
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is proved using variational inequalities [7]. The case when the quadrature 



formula (1.2) involves derivatives but is supported at only one point is espe- 
cially relevant to the present paper and is discussed in Section 3. In each of 
these cases, the existence theorems are true in any number of dimensions. 

In the plane, quadrature domains can be constructed explicitly. The only 
explicit examples in higher dimensions are a sphere in IR n and some special 
examples in M 4 discovered by L. Karp [12]. 

The current investigation is motivated by the lack of explicit examples in 
the higher-dimensional case and the lack of qualitative understanding. We 
address the question Are quadrature domains algebraic in M. n , with n > 3? 
which was raised by H. S. Shapiro [TO] (p. 40) and posed again as an open 
problem by B. Gustafsson [I] (p. xv). 

As usual, by "algebraic domain" we mean that the boundary of the do- 
main is contained in the zero-set of a polynomial. 

In the plane, under no regularity assumptions on dQ, it was shown by 
D. Aharanov and H. S. Shapiro (1976) that quadrature domains are always 
algebraic. B. Gustafsson (1983) showed that they have further nice properties 
[8] which we summarize in the following. 

Theorem 1.1 (B. Gustafsson, 1983). If Q C 1R 2 is a quadrature domain, 
then Q is algebraic. Moreover, for some polynomial P(x, y), the boundary dQ 
consists of all of the points of{P(x,y) = 0} except finitely many. Moreover, 
the leading order (homogeneous) term in P is a constant times (X 2 +Y 2 ) M for 
some M. Moreover, there are some explicit relations between the coefficients 
of P and the coefficients a m ^ in the quadrature identity. 

By the remark above, this applies as well to quadrature domains for 
harmonic functions in the plane, so it makes sense to ask if any part of the 
theorem extends to M n . Considering that the examples in M 4 constructed 
by L. Karp [12] were algebraic, one might hope that the answer is "yes". 
We show that in fact quadrature domains are not always algebraic, thus 
confirming the following Conjecture [TS] in n = 4 dimensions. 

Conjecture 1.2. In all dimensions greater than two, there exist quadrature 
domains that are not algebraic. 

Our answer is based on constructing explicit examples in Section |3j In 
Section [4] we use these same examples to generate some exact solutions to 
the Laplacian growth problem. 

First, we review the definition of the Schwarz Potential in the next section. 
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2 The Schwarz Potential 



Suppose that T is a non-singular, real-analytic curve in the plane. Then the 
Schwarz function S(() is the function that is complex-analytic in a neigh- 
borhood of T and coincides with ( on T (see |2] for a full exposition). If T 
is given algebraically as the zero set of a polynomial P(x,y), by the implicit 
function Theorem we can obtain S(Q by making the complex-linear change 
of variables ( = x + iy, ( = x — iy, and then solving for ( in the equation 

P I — - — , —z— ) = (See [HI p. 3]). On the other hand, S(() can also be 

written in terms of the conformal map / from the unit disk to the domain, 
using the formula 

s(0 = r(i^), (2.i) 



where /* denotes the function obtained by conjugating the coefficients of /. 

In the next section, we will utilize both such representations of S(Q, so 
let us illustrate each by an example. 

Example ("C. Neumann's oval"): Let a > be a real parameter. Sup- 
pose that T is the Jordan curve defined by 

(x 2 + y 2 ) 2 = a 2 (x 2 + y 2 ) +4x 2 . 

See Figure ([TJ, for a plot of T for different values of a > 0. 




Figure 1: A plot of T for different values of a > 0. 
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Changing the variables we have (CC) 2 = fl2 (CC) + (C + C) 2 - Solving for ( 
gives 

q(n = C(a 2 + 2) + 2Cy/a 4 + a 2 + C 2 
lU 2(C 2 -1) 

indicating that T is the boundary of a quadrature domain since S(() is mero- 
morphic in the interior with only poles at ( = ±1. 



Let us make the same observation using the formula (2.1 ). The conformal 
map from D to the interior of T is given by 

/«- {Rt - 1)z 



R(R 2 - z 2 ) ' 



a _|_ ^J a 2 _)_ 4 

where i? = > 1. Without actually calculating the inverse of /, 

but just observing that it is a conformal map from the interior of V to D, 

we can see that S{() = f* J * s meromor phic in the interior of T and 

has two simple poles at f(±l/R) = ±1. 

Suppose that T is more generally a nonsingular analytic hypersurface in 
W 1 , and consider the following Cauchy problem posed near T. The solution 
exists in some neighborhood of T and is unique by the Cauchy-Kovalevskaya 
Theorem. 

Aw = near V 

Hr = ^||x|| 2 (2.2) 
Vw\r = x 



Definition. The solution u>(x) of Cauchy problem 2.2 is called the Schwarz 
Potential of T. 

In IR 2 , the Schwarz function can be directly recovered from the Schwarz 
potential. Consider S(Q = 2d^w = w x — iw y . The Cauchy-Riemann equa- 
tions for S follow from harmonicity of w, and Vw = x on T implies S(() = ( 

on r. 

As mentioned in the introduction, quadrature domains have an equiva- 
lent definition framed in terms of the Schwarz potential. Namely, the follow- 
ing theorem shows that condition (ii) in the Introduction has a counterpart 
in higher dimensions, where instead of "meromorphic", we have that the 
Schwarz potential has a single-valued real-analytic continuation throughout 
Q except for finitely many points, where Aw is a distribution of finite order. 
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Theorem 2.1 (Khavinson, Shapiro, 1989). Suppose that Q is a connected 
domain in M. n with C 1 boundary. Then Q is a quadrature domain with quadra- 



ture formula (1.2) if and only if the Schwarz potential of dVL satisfies 

w(ar)=n(ff(x)-Q(x)), (2.3) 
where H(x) is harmonic in Q and C 1 in Q, and Q(x) is the result of applying 



12-71 



the quadrature formula l[1.2\) to the fundamental solution E n (x.,y) = c n ||x — 

jx|| 2 = r 2 } be a sphere of radius r. When 



(in its second argument). 



Example: Let T := {x G I 

n — 2, it is easy to verify that 



w(()=r 2 (log|C|+l/2-log(r)) 



solves the Cauchy Problem (2.2), and in higher dimensions the Schwarz po- 
tential is 



n 



w x 



-r 2 - 



2(n-2)' (n-2)||x|| n - 2 ' 

Let us check that this agrees with Theorem |2.1| The mean- value property 
for the ball of radius r about zero, -B r (0), gives the quadrature formula: 

udV = Vol (B r )u(0). 

B r (0) 

Applying this quadrature formula to E n (x, y ) in its second argument gives 
Q(x) = Vol(5 r )£ n (x) = Vol^^CnllxH 2 -" = ^VoKSOcnllxH 2 -" 

„n 

1 1 1 12 — « 



n(n — 2) 

since the constant c n appearing in the fundamental solution is 
where u n is the surface area of the unit sphere. Taking H(x) 



(n - 2)ujn 
n 



2 n - 2) 



we obtain the same formula (2.3). 
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3 An example of a non-algebraic quadrature 
domain 



A simple case of particular interest is when the quadrature formula (1.2) is 
supported at a single point. 



M 

udV = a a d a u(0), (3.1) 
M=o 

where a = (cui, 02, ..,«„) again is a multi- index. 

This quadrature formula corresponds to a domain that has only finitely 
many non- vanishing harmonic moments. In two dimensions, this is equivalent 
to the conformal map from the disk being a polynomial. 

Under the condition that the coefficients a a for \a\ > are small in 
comparison with the leading coefficient a(o,o,..,o)> the existence and uniqueness 



of domains admitting a quadrature formula of the type (3.1 ) follows from one 



of the earliest results on the solvability of the inverse potential problem [20J 

(cf. M). 

Let us further specialize to the case when all derivatives in the quadrature 



formula (3.1) are with respect to the same variable, say x\. 



udV = $>fc^(0). 

fc=0 1 



(3.2) 



In [12], special cases were explicitly generated by rotating the limacon 
D :={(:( = w + aw 2 , \w\ < 1} about the x-axis, where a > 1/2 is a real 
parameter. Before rotation, this starts out as a two-dimensional quadrature 
domain with quadrature formula 



f f (()dA = tt[(1 + 2a 2 ) f(0) + a f 
Jd 



The four-dimensional domain generated by rotation is a quadrature do- 
main with quadrature formula 



udV 



7l~ 



n 



1 + 6(x 2 + 2a 4 



2a 2 ) du , . a 2 d 2 u , ' 
■(0) + 77^(0) 



dxi 



\2dx\ 



where x\ corresponds to the axis of symmetry. 

£ 4 by constructing a domain admitting the 



We prove Conjecture 1.2 



in 



quadrature formula (3.2) with m = 1, ciq > 0, and a\ > 0. 
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Theorem 3.1. In M 4 , there exist quadrature domains that are not algebraic, 
namely, admitting the quadrature formula (3.2) with m — 1, do > 0, and 
cl\ > (with a\ small in comparison with oq). Moreover, the boundary can 
be described explicitly in terms of elliptic integrals of the third kind. 

We will work in terms of the Schwarz potential introduced in the previous 

section. Before beginning the proof we make an observation regarding axially 

symmetric potentials. 

Recall the axially symmetric reduction of Laplace's equation: 

Suppose that u(x±, x 2 , x n ) is harmonic in M. n and axially symmetric 

about the xi-axis. Write U(x,y) = u(xi, x 2 , x n ), where x = x\ and y = 

\fx1 + ■■■ + x 2 n . Then, 

y 

Indeed, Am = div (U x , U y y X2 ,U y y X3 , U y y Xn ) 



u xx + Uyy[J2yli) +u y[J2 y*** 



i=2 I \i=2 



where easy calculations give Y^=2 yli = ^ an< i ^ x » a! < 



(n-2) 



y 

2U 

In the case n = 4, U satisfies the equation AU H = 0, iff yU (x, y) is 

y 

a harmonic function of two variables. Indeed, 

A(yU) = yAU + 2W ■ Vy + UAy = yAU + 2U y . 

This fact about axially symmetric potentials in IR 4 was used in [13] and [12] 
where it is explained in more detail. 

If Q is a domain in M 4 with axially symmetry, then the above consider- 
ations apply to the Schwarz potential, since the rotational-symmetry of the 
Cauchy data (2.2) is passed to the solution w. 

Thus, the Cauchy problem (2.2) defining the Schwarz potential w is 
reduced to the following two-dimensional Cauchy problem for W(x,y) = 
w(x 1 ,x 2 ,x 3 ,x i ): 

AW + ^—^ = 0, near T 

y 

W| 7 = (x,y) 
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where 7 is the symmetric curve in the plane whose rotation generates the 
boundary of fl. 

According to the previous observation, we have 

A(yW(x,y)) = 0, (3.4) 

which will be used below. 



Proof of Theorem 3. 1 , We will construct fl by first describing a conformal 
map / : D — > D p from the unit disk D to a domain D p in the plane symmetric 
with respect to the real axis. Then we will take fl to be the domain generated 
by rotation of D p into M 4 about the x-axis. i.e., (x\, 22, £3, £4) £ if and only 

if ^xi, \/ x\ + x\ + x"lj G D p . Notice that the boundary T = dQ is algebraic 

if and only if 7 = dD p is algebraic, in which case 7 = {p(x, y) = 0}, and 

r = ip (xi, \/ x\ + x\ + x\ J = o|, where p(x, y) is a polynomial that is even 

in the variable y. 

The conformal map / and its relevant properties will be established in 
the proofs of the following claims. 

Claim 1: There exists a function / with real coefficients, analytic and uni- 
valent in a neighborhood of D, with /(0) = 0, so that the function 

=0(*) (3-5) 

is rational, and its only pole in D is a pole of exact order 3 at z = 0. 
Claim 2: The (analytically continued) function / has an infinitely-sheeted 
Riemann surface. 

We defer the proofs of the claims in favor of first seeing how they are used 
to prove the theorem. 

Since, by Claim 1, D — /(D) is the image of D under a univalent function 
analytic in a neighborhood of D, the boundary of D p , and therefore of fl, is 
analytic and non-singular. Thus dfl has a Schwarz potential w(xi, x 2 , x 3 , x 4 ), 
which by axial symmetry can be reduced to a function W(x,y) of two vari- 



ables (see the discussion above just before the proof). As stated in (3.4), 
yW(x,y) is harmonic near dD p . 

Let ( = x + iy. Then the function 

V(x,y)=yW(x,y)-%{( 3 }/12, 
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is harmonic (subtracting the harmonic term ^s{( 3 }/12 simplifies the following 
calculations without changing the singularities of V). 

Next, consider the (complex-analytic) function obtained by taking the 
(■-derivative = \{V X — iV y ). 



The Cauchy data for W, stated in (3.3), can be written W = CC/2 and 
W( = C/2, so that on the boundary 7 = dD p , the equation above gives 

V\ + _ * ( t A 2 

v ch--^ + -^\Q + -^-i\s-Q > 
C-C 

where we have also used y = . 

Now we replace ( with S(() the Schwarz function of 7. Then, both sides 
of the equation are analytic in a neighborhood of 7, and thus the following 
becomes an identity (not just valid on the boundary). 

^C = ^(C-S(C)) 2 . (3.6) 

This relates the 4-dimensional Schwarz potential to the 2-dimensional 
Schwarz function. 

As follows from (2.1) when the coefficients of / are real, the Schwarz 
function of 7 satisfies the functional equation 

*(/(.-.)) = /: 1 



Using this relationship, and substituting ( = f{z) into (3.6), we get 



■ / / 1 \ \ 2 



V C (f(z)) = l [f(z)-f 



By Claim 1 the right hand side is g(z), a function analytic in a neighborhood 
of D, except for a pole at z = of order 3. So, g(z) = — with Q(z) 
analytic in a neighborhood of D and Q(0) 7^ 0. 
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Since /(O) = and / is univalent, we have / = CMO> ^ or some MO 
analytic and non-vanishing in D p . Thus, 

vdO - g[t (0) - , 

which is analytic in D p except for a pole of order 3 at £ = 0. 
This implies that (see [121 P- 183]) 

w(xi,x 2 ,x 3 ,x 4 ) = iu(x) = A ■ ||x||~ 2 + 5 ■ 9 Xl (||x|| -2 ) + #(x), 

where A and i? are constants and H(pc) is harmonic. 

Since, up to a constant, ||x|| -2 is the fundamental solution Ki(x, 0) evalu- 
ated at y = 0, Theorem 2.1 implies Q is a quadrature domain with quadrature 
formula of the form 



r 0u 

/ udV = a u{0) + ai — (0). 



Next we apply Claim 2 in order to show that the boundary of Q is not 
algebraic. 

Note: It does not follow immediately from Claim 2 that the image of <9B 
is non-algebraic, since it is possible for an algebraic domain to be conformally 
mapped from the unit disk by a transcendental function. Indeed, any non- 
circular ellipse is such an example. 

Suppose, to the contrary, that 7 = dD p is the zero set of a polynomial. 
Then, as observed in Section |2j the Schwarz function S(Q is an algebraic 
function, since it can be obtained by solving for ( in the equation 

P ( ^> , I 0. 

2 ' 2i 



Substitute z = f into (3.5) and use the formula (2.1) 



/-HO 



where again, the conjugation is missing because the coefficients of / are real. 
We obtain 

g(r l (0) = (C-S(C)) 2 - 
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Since g(z) is rational, and S(Q is algebraic, this implies 

r 1 (0 = g- 1 (((-s(()) 2 ) 

is also an algebraic function. This contradicts the fact that / _1 (C) is the in- 
verse of a transcendental function. We conclude that S(() is transcendental, 
and the boundary of Q is not contained in the zero set of a polynomial. 

□ 

It remains to prove the claims. 
Proof of Claim 1. With a > 0, C > real parameters, let 



h(z) = C^Ap^l + az). 
z V z 

It is easy to see that the square root above has two single valued analytic 
branches on \z\ = 1. We choose that branch which is > at z = 1. 
We observe that 

h(l/z) = -h{z). (3.7) 

Since h(z) is analytic in a neighborhood of <9B, say {1 — e < |z| < 1 + e} 
it has a Laurent expansion 



Kz) = 



Cj z 



j 



Let f(z) = Y^jLi c j z ' By (3.7), Cj = —Cj, and c = 0. 
Therefore, 



h(z) = f(z) -/(-), (3.8) 



where f(z) is analytic in a neighborhood of D, and f(l/z) is analytic in a 
neighborhood of C\B with f(l/z) — > as 2; — > 00. 
Observe that /(z) satisfies (3.5) if we take 



2° 



We remark that, up to changing the sign of the parameter a, this is the only 
choice of g for which all the conditions in the claim can be satisfied. 



12 



For toeD, f(w) can be expressed as the contour integral 



2iri J\ z \=i z — w 



Indeed, by (3.8) 



\ Z \ =1 Z-W J ]z \ =1 Z-W J ]z]=1 Z-W J ]z \ =1 Z-W 

and we have / dz = 0, since the integrand is analytic in C\ D and 

J\ z \=i z-w 

is 0(l/|z| 2 ) as z — > oo. Moreover, since f(z) is analytic in a neighborhood 

of D, the Cauchy integral formula gives / dz = f(w). 

2m J ]zl=1 z-w 



So we have 



f( w ) = / - ( S dz - ( 3 - 9 ) 

2ni J\ z \ =l z(z-w) 

Next, we show that f(w) is univalent at least for sufficiently small values 
of the parameter < a < 1. The parameter C does not affect the univalence 
of f(w), so we fix C = 1. 

Use (3.9) to write the derivative f'(w): 

/» = - [ + + 

2wi J\- A=1+£ z(z - w) 2 

where we have also deformed the contour slightly with e > 0. As a — > 0, the 

- 1) 

integrand converges (uniformly for \z\ — 1 + e) to — -r. Thus, f'(w) 



converges uniformly in D to 
1 



z(z — w)< 



-dz 



2iri J\ g \ =1+e (z - w) 2 z(z - w) 
using residues. 

This guarantees univalence of f{w) throughout D for all sufficiently small 
a > 0. Indeed, for any two points w± and w^ in D, 

f(w 2 ) - f(wi) = / f'(w)dw, 
13 



and f'(w) is uniformly close to 1 for a > sufficiently small. 

□ 

Proof of Claim 2. We investigate the global analytic function whose princi- 



pal branch is given by the integral (3.9) 



We will deform the contour in (3.9), but first we manipulate the integrand 



so that the new integral will converge. Using 



[z 2 — 1) w w 2 



w 



z(z — w) z z(z — w) 

we have 

f(w) = A + A x w + C(w 2 - l)F{ 
where A Q and A\ are integrals that do not depend on w, and 

p(w)= if V(i + f)(i + ^ 



2m J ]z]=1 z(z - w) 

Now deform the contour until it it "just surrounds" the segment (— oo, — 1/a]. 
The integrand changes sign as it switches sides, and the path changes direc- 
tion. Thus, we obtain twice the integration along a segment: 

m y_ 00 z{z - W) TTl (z - w) 

where G(z) = (l+-)(-+a)-. By further manipulation, this can be expressed 
in terms of a standard form of a complete elliptic integral of the third kind, 
which we will do at the end of this section after the proof. For now, we show 



directly from (3.11) that the analytic continuation of F(w), and hence f(w), 
has an infinitely- sheeted Riemann surface. 

The monodromy for F(w) can be obtained as a nice exercise in applying 
the Sokhotski-Plemelj relation [TUJ Ch. 14]. Accordingly, the jump across 



the segment (— oo, — 1/a] of such an integral as appears in (3.11) is equal to 
1i\i times the numerator of the integrand, i.e., let w approach a point x on 
the segment (— oo, —1/a] from above and below, respectively, then 



2^fG(x). (3.12) 

In order to perform analytic continuation, choose a base point Wq, say 
in the upper half-plane, and let Fq(w) denote the principal branch (which is 
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expressed by the integral itself). Then, let 71 be a loop based at w that winds 



once around I /a. According to (3.12), performing analytic continuation along 
71, we obtain a new branch Fi(w) = F (w) + 2y/G(w). If 71 is traced again, 
then we return to the original branch since 2yjG(w) switches sign and cancels 
the contribution from Fq. So w = 1/a only has ramification index 2. 

Consider a second loop 7 2 based at w that winds around the origin, 
crossing the segment [0, a]. Traveling along this loop does nothing to the 
initial branch Fq, but if we first continue along 71, then the branch F% will 
be affected when we continue along 72. Alternate between these two loops 
and perform analytic continuation along 71 • 72 • 71 • • • 71, where 71 is traced 
k + 1 times. Let F^w) denote the new branch obtained at each step. Then 
we have 

F 1 (w) = F (w) + 2 v fGjzj 

F 2 (w)=F (w)-2^/Gjz) 
F 3 (w) = F Q (w)+2^Gjz) + 2 v ^Gjz) 

F 4 (w) = F (w) - 
F 5 (w) = F (w) + 2 ^Gjz) + A^Gjz) 
F 2k+1 {w) = F (w) + 2(k + l)y/G{z) 
Thus, the Riemann surface for F(w) is infinitely-sheeted. 

□ 

Reduction of the integral to a standard form: Let us make f(w) 
more explicit by red ucing F(w) to a standard form. Make the substitution 



■i dz=§, in (3.11 



to get 



Using —j - j^-, 



7T 7 (c+l» 

where Cq is an integral that doesn't depend on w. 
For the remaining integral we write 
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ra y/(t-l/a)(t- a )t = ( j )( _ } ra 1 d j. 

JO (e+l/w) s / AS Jo (^+l/-!«) A /(5-l/a)(e-a)5 S 

This integral on the right hand side is entry 4 from section 3.137 of [6], 
where it is expressed in terms of a complete elliptic integral of the third kind. 
The argument w appears only within the so-called elliptic characteristic. 



4 Laplacian growth in 



Given an initial domain Q Q , consider the following moving-boundary problem. 
Find a one-parameter family of domains, {fit} C M. n so that the normal 
velocity, v n , of the boundary T t := dQ t is determined by Green's function, 
P(x, t), of Qt with a fixed singularity positioned at xo G Qf 

v n \r t = -VP 
— n ;™ o. 

(4.1) 



AP = 0, in Q t 



P\r t = 

P(X X , t) ~ Q ■ P n (x, X ), 

where E n is the fundamental solution of the Laplace equation with singularity 
at xo and Q > (Q < 0) determines the injection (suction) rate at the source 
(sink) x . 

If the domains fl t are bounded as in the case considered below, with 



Q > problem (4.1) actually produces a shrinking boundary. We get a 
growth process if fl t contains infinity so P then solves an exterior Dirichlet 
problem. In such a situation, it is common to place the sink at infinity 
by prescribing asymptotics for VP so that the flux across neighborhoods of 
infinity is proportional to Q. 

This is a non-linear moving boundary problem that models viscous fin- 
gering in a Hele-Shaw cell (when n = 2) or bubble growth in a porous media 
(when n = 3). These processes exhibit complicated pattern formation |16j . 
yet miraculously, there turns out to be an abundance of explicit, exact solu- 
tions in the plane. This miracle is partly explained by Richardson's Theorem 
[IB], which guarantees that the property of fl being a quadrature domain is 
preserved under Laplacian growth, and moreover the time-evolution of the 
quadrature formula is simple. 



Theorem 4.1 (S. Richardson). If Vt t solves Problem (4-1), then for any 
harmonic function u 

,1 r 

udV = Qu(x ). 



d_ 
dt 
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Given the many equivalent definitions of quadrature domains a few of 
which were mentioned in the introduction, it is not surprising that there are 
alternative statements of Richardson's Theorem. For instance, applying the 
Theorem to a basis of harmonic polynomials expanded about the point xq, we 
obtain that all harmonic moments are constant except one (infinitely-many 
conservation laws). One may ask if there is a moment-generating function for 
the harmonic moments. The answer is that it is precisely the exterior grav- 
itational potential of fit, which according to Richardson's Theorem evolves 
by inheriting the effect due to a single additional point-mass at xo growing 
linearly in strength. 

Yet another interpretation is relevant to the approach taken in the above 
sections. Namely, all singularities of the Schwarz potential are stationary ex- 
cept one positioned at xq which does not move but simply grows in strength. 



This is a consequence of Richardson's Theorem combined with Theorem 2.1 



It can be seen more directly by establishing the following formula, shown in 
[T5] from elementary calculations, relating the Schwarz potential w to the 
pressure P. 

^w(x,t) = -nP(x,t), (4.2) 

where n is the spatial dimension. 

In the statement of the problem (4.1), if we allow P to have more exotic 
singularities, then the quadrature domains constructed in [T2] provide exact 
solutions in M 4 with a dipole flow superimposed on a source of varying rate. 
Some exact solutions without multipole flows were described in [15] , but they 
required a combination of sources and sinks. The examples constructed in 
Section [3] provide exact solutions to the Problem (4.1) as stated (i.e. with 
a single sink). Before describing the solution, we review its analogue in the 
plane, which is the well-known example (due to Polubarinova-Kochina |17j ) 
that encounters a cardioid. 

Example 1 (IR 2 ): Consider the family of domains D with boundary given 
by the curves {( : ( = az 2 +bz, \w\ < 1} with a > 0, b > 2a real. The Schwarz 
function is given by S(() = -2ab/{a - a/o 2 + 4b() + 4b 3 /(a - ^ a 2 + 4&C) 2 
which has a single-valued branch in the interior of the curve for appropriate 
parameter values a and b. The only singularities of the Schwarz function inte- 
rior to the curve are a simple po le an d a pole of order two at the origin. Since 
S(Q = 2d c w(x,y,t), Equation becomes |5(C,t) = -Ad c P(x,y,t). So, 



in order for a one-parameter family of domains to solve Problem (4.1), the 
singularities of S(() must be time-independent except for one simple pole. 
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Given an initial domain from this family we can choose a one-parameter slice 
of domains so that the simple pole increases (resp. decreases) while the pole 
of order two does not change. This gives an exact solution to the Laplacian 
growth problem with injection (resp. suction) taking place at the origin. In 
the case of injection, the domain approaches a circle. In the case of suction, 
the domain develops a cusp in finite time. 

Example 2 (IR 4 ): Similarly, the the domains constructed in the previ- 
ous section contain one-parameter families that solve the Laplacian growth 
problem with a sink at the origin. Since changing C simply scales the do- 
main, the different shapes that can occur in these one-parameter families are 
determined by the value of a. Near the value a ~ 0.82217..., the boundary 
develops a cusp (apparent in Fig. [2]). 




Figure 2: Two-dimensional profile of a solution of Laplacian growth in M 4 , 
reminiscent of the cardioid example in the plane. 



5 Concluding remarks 

1. Before rotation into M 4 , the domain D p is not a quadrature domain, but 
it is a quadrature domain in the wide sense and admits a formula for the 
integral of any function F analytic in D p in terms of a distribution supported 
on the interval [— a, 0], where —a = /(—a) is the image of —a under the 



conformal map described in the proof of Theorem 3.1 
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2. It is possible to obtain an explicit (but complicated) formula for the 
Schwarz potential of an axially symmetric domain in IR n . First, take the 
axially symmetric reduction to (two variables) as discussed just before the 



proof of Theorem 3.1 



. (n - 2)U X 
AU+- — = 0, 

x 

where U (x, y) is the Schwarz potential in the cylindrical variables, and we 
have taken y as the variable along the axis of symmetry in order to correspond 
with a reference used below. 

Now, let x and y each take complex values and make the change of vari- 
ables to characteristic coordinates z = x + iy , w = x — iy: 

(n-2)(U z + U w ) _ 
Uzw+ 2(z + w) " U - 

This is a complexified hyperbolic equation with lower order terms. The 
Cauchy problem can be solved using a generalization of d'Alembert's formula 
involving a so-called Riemann function A(s,t;z,w). The general procedure 
is discussed in (5J Ch. 4 and Ch. 5], and the Riemann function for this 
equation is given by formula (5.36). 

A(s , t ., z , w) = . . (' + «>' r (a i, 1; (• - ^ - »> 



(s + w) x / 2 (t + z) x / 2 V 2 ' 2 ' 1 (s + w)(t + z) 

where F(a, b; c; x) denotes the Gauss hypergeometric function, and A = (n — 
2). In the case n — 4, the hypergeometric function F(l,l;l;;r) = is 
rational, but in the case n = 3 it is transcendental. 

Let S(z) be the Schwarz function of the curve 7 that generates the surface 
(by rotation). The complexification of the curve then becomes the graph in 
C 2 of S(z). Using also the fact that the inverse S~ 1 (w) = S(w), is just the 
function obtained by conjugating the coefficients of S, the following formula 
for U(z, w) is a complex version of the formula (4.73) from [5]: 

4U(z, w) = 2 [A(S(w), wj z } w)S(w)w + A(z, S(z); z, w)_S(z)z] 
+ fs (z) ^r t A(S{t),t; z, w)S{t)t+A(S(t),t; z, w)S(t)-A t (S(t),t; z, w)S(t)tdt 

-l! {w) 1 ^A(a,S{s);z,w)S(s)a 

+A(s, S(s); z, w)S(s) — A s (s, S(s); z, w)S(s)sds, 
where we have also taken into account the Cauchy data in the above formula: 

" £/| 7 = zw/2 
U z \~ ( = w/2 

U w \ry Z /2 
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Instead of using the Schwarz function, this formula can be written in terms 
of the conformal map / from the disk by making the change of variables t = 
s = f(v)^ an d using the relation S(f(rj)) = /Q)- Making this change of 
variables, one can then consider the following problem whose solution would 
yield three-dimensional quadrature domains with prescribed quadrature: 

Suppose the singularities of U(z,z) are prescribed throughout Q. Is it 
possible to then determine / using the formula above? 

Conceptually, this is analogous to what was done in Section [3j recovering 



f(z) from the functional equation (3.5), but in practical terms it appears 
much more difficult. 

3. Since M 4 has proven easier to work with than the more relevant M 3 , one 
wonders if axially symmetric examples in IR 3 can be expected to "interpolate" , 
in some sense, between ones in R 2 and corresponding ones in R 4 , so that one 
can formulate conjectures regarding the more physically relevant case of M 3 . 
For instance, starting with a circle, a sphere in M 3 , and a hypersphere in 
R 4 (all of the same radius), consider a dipole flow from the center (in the 
direction of the ) . The two and four-dimensional examples each develop 
a cusp on the x-axis. We expect the three-dimensional example does also, and 
it is tempting to conjecture that the position of the cusp is located between 
the positions of the two and four- dimensional cases. It may be more suitable 
to formulate this question generally within the framework of "generalized 
axially symmetric potentials" [21], i.e. considering non-integer values of the 
parameter A in the equation 

Au + — = 0. 

x 

Then the question is whether certain aspects of Laplacian growth "depend 
monotonically" on A. 



4. The procedure used in the proof of Theorem 3^ can also be used to 
construct other domains having quadrature formula of the form (3.2) with 
any choice of m, i.e. having any number of non-vanishing harmonic mo- 
ments. Considering these other cases, one observes that only exceptional 
cases are algebraic. It is tempting to conjecture that quadrature domains 
are almost never algebraic in dimensions greater than two, and perhaps this 
would help explain why it has been so difficult to obtain exact solutions to the 
higher-dimensional Laplacian growth problem. For Laplacian growth in any 
number of dimensions, if the domain is initially a quadrature domain, then 
as mentioned in Section [4] it remains a quadrature domain, but in the plane 
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it moreover remains an algebraic curve of bounded (uniformly throughout 
time) degree. 

5. Instead of prescribing the coefficients in the quadrature formula and 
trying to find the domain, one can consider a given domain, say with algebraic 
boundary, and try to determine if it admits a quadrature formula. In this 
case, it is more reasonable to consider "quadrature domains in the wide sense" 
and allow for a measure supported on a continuum instead of a finite sum of 
point-evaluations. For an analytic surface, it is always possible to find such 
a measure supported on a slightly smaller domain. For ellipsoids, there is a 
measure supported on a two dimensional set inside called the "focal ellipse" 

Acknowledgement: We wish to thank Lavi Karp for directing us to an 
important reference |20j . 

References 

[1] D. Aharonov, H. S. Shapiro, Domains on which analytic functions satisfy 
quadrature identities, J. Anal. Math. 30 (1976) 39-73. 

[2] P. J. Davis, The Schwarz function and its applications, The Mathematical 
Association of America, Buffalo, N. Y., 1974. The Carus Mathematical 
Monographs, No. 17. 

[3] B. A. Dubrovin, S. P. Novikov, A. T. Fomenko, Modern geometry: meth- 
ods and applications, part I, Moscow 1986. 

[4] P. Ebenfelt, B. Gustafsson, D. Khavinson, M. Putinar (eds.), Quadrature 
domains and their applications, Oper. Theory Adv. Appl., 156, Birkh 
user, Basel, 2005. 

[5] P. Garabedian, Partial differential equations, 2nd ed., Chelsea Pub. Co., 
1998. 

[6] I. S. Gradshteyn, I. M. Ryzhik, Table of integrals, series, and products, 
Academic Press (English translation edited by A. Jeffrey), 1965 

[7] B. Gustafsson, Applications of variational inequalities to a moving bound- 
ary problem for Hele-Shaw flows, research report TRITA-MAT-1981-9, 
Royal Inst, of Technology, Stockholm. 



21 



[8] B. Gustafsson, Quadrature identities and the Schottky double, Acta Ap- 
plicandae Math. 1 (1983), 209-240. 

[9] B. Gustafsson, A distortion theorem for quadrature domains for harmonic 
functions, J. Math. Anal. Appl. 202 (1996), 169 - 182. 

[10] P. Henrici, Applied and Computational complex analysis, vol III, John 
Wiley and Sons, 1986. 

[11] V.K. Ivanov, An inverse problem for the potential of a body close to a 
given one, Izvestiya AN SSSR, Ser: Math. 20 (1956), 793-818 (Russian). 

[12] L. Karp, Construction of quadrature domains in W 1 from quadrature 
domains in R 2 , Complex Var. Elliptic Eq., 17 (1992), 179-188. 

[13] D. Khavinson, On reflection of harmonic functions in surfaces of revo- 
lution, Complex Variables 17 (1991), 7-14. 

[14] D. Khavinson, H. S. Shapiro, The Schwarz potential in W 1 and Cauchy's 
problem for the Laplace equation, TRITA-MAT-1989-36, Royal Institute 
of Technology, Stockholm. 

[15] E. Lundberg, Laplacian growth, elliptic growth, and singularities of the 
Schwarz potential, J. Phys. A: Math. Theor., 44 (2011), 135202. 

[16] J. Mathiesen, I. Procaccia, H. L. Swinney, M. Thrasher The Universality 
Class of Diffusion Limited Aggregation and Viscous Fingering, Europhys. 
Lett., 76 (2006), 257-263. 

[17] P. Ya. Polubarinova-Kochina, On the motion of the oil contour, Dokl. 
Akad. Nauk. S.S.S.R. 47 (1945), 254-257 (in Russian). 

[18] S. Richardson, Hele-Shaw flows with a free boundary produced by the 
injection of fluid into a narrow channel, J. Fluid Mech., 56 (1972), 609- 
618. 

[19] H. S. Shapiro, The Schwarz function and its generalizations to higher 
dimensions, Wiley-Interscience, 1992. 

[20] L. N. Sretenskii, On the inverse problem of potential theory, Izvestiya 
Acad. Nauk SSSR, Ser. Mat. 2 (1938), 551-570 (Russian). 



22 



[21] A. Weinstein, Generalized axially symmetric potential theory, Bull. 
Amer. Math. Soc. 59 (1953), 20-38. 



Purdue University 

West Lafayette, IN 47907 

eremenko@math.purdue. edu 

elundber@math.purdue.edu 



23 



